(* Content-type: application/mathematica *)

(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)

(* CreatedBy='Mathematica 6.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       145,          7]
NotebookDataLength[     71165,       1922]
NotebookOptionsPosition[     65285,       1722]
NotebookOutlinePosition[     65657,       1738]
CellTagsIndexPosition[     65614,       1735]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{
Cell[BoxData[{
 RowBox[{
  RowBox[{"Off", "[", 
   RowBox[{"ClebschGordan", "::", "\"\<tri\>\""}], "]"}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Off", "[", 
   RowBox[{"ClebschGordan", "::", "\"\<phy\>\""}], "]"}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Off", "[", 
   RowBox[{"SixJSymbol", "::", "\"\<tri\>\""}], "]"}], ";"}]}], "Input",
 CellChangeTimes->{{3.456663504453125*^9, 3.45666350465625*^9}, 
   3.54004256367356*^9}],

Cell[CellGroupData[{

Cell["Physical Constants", "Subtitle",
 FontColor->RGBColor[0, 0.356863, 0]],

Cell[BoxData[{
 RowBox[{
  RowBox[{
  "oFile", "=", 
   "\"\<C:\\\\Users\\\\work\\\\Data\\\\output\\\\Rb87-D1.txt\>\""}], ";", " ", 
  RowBox[{"(*", 
   RowBox[{"output", " ", "file"}], "*)"}], "\[IndentingNewLine]", 
  RowBox[{
  "pathname", "=", "\"\<C:\\\\Users\\\\work\\\\Data\\\\2012-05-15\\\\\>\""}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"filename", "=", "\"\<Rb-D1-line.dat\>\""}], ";", " ", 
  RowBox[{"(*", "datafile", "*)"}], "\[IndentingNewLine]", 
  RowBox[{"ExpSpecFile", "=", 
   RowBox[{"pathname", " ", "<>", " ", "filename"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"dataTab", "=", " ", 
    RowBox[{"Import", "[", 
     RowBox[{"ExpSpecFile", ",", "\"\<Table\>\""}], "]"}]}], ";"}], 
  "\[IndentingNewLine]", 
  RowBox[{"(*", 
   RowBox[{
   "Form", " ", "repetition", " ", "rate", " ", "list", " ", "and", " ", 
    "offset", " ", "frequency", " ", 
    RowBox[{"list", ".", " ", "Extract"}], " ", "from", " ", "actual", " ", 
    "data", " ", "file"}], "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"fRepList", "=", 
   RowBox[{"{", "}"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"fOffsetList", "=", 
   RowBox[{"{", "}"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"n", "=", "1"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"While", "[", 
   RowBox[{
    RowBox[{"n", "\[LessEqual]", " ", 
     RowBox[{"Length", "[", "dataTab", "]"}]}], ",", "\[IndentingNewLine]", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"fRepList", "=", 
       RowBox[{"Append", "[", 
        RowBox[{"fRepList", ",", 
         RowBox[{
          RowBox[{"dataTab", "[", 
           RowBox[{"[", "n", "]"}], "]"}], "[", 
          RowBox[{"[", "2", "]"}], "]"}]}], "]"}]}], ";", 
      RowBox[{"(*", 
       RowBox[{"\"\<2\>\"", " ", "=", " ", 
        RowBox[{
        "2", "nd", " ", "column", " ", "of", " ", "the", " ", "datafile", " ",
          "contains", " ", "the", " ", "rep", " ", "rate"}]}], "*)"}], 
      "\[IndentingNewLine]", 
      RowBox[{"fOffsetList", "=", 
       RowBox[{"Append", "[", 
        RowBox[{"fOffsetList", ",", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"-", "2"}], "*", 
           RowBox[{
            RowBox[{"dataTab", "[", 
             RowBox[{"[", "n", "]"}], "]"}], "[", 
            RowBox[{"[", "1", "]"}], "]"}]}], ")"}]}], "]"}]}], ";", " ", 
      RowBox[{"(*", 
       RowBox[{"\"\<1\>\"", " ", "=", " ", 
        RowBox[{
        "1", "st", " ", "column", " ", "of", " ", "the", " ", "data", " ", 
         "file"}]}], "*)"}], "\[IndentingNewLine]", 
      RowBox[{"n", "++"}], ";"}], "}"}]}], "]"}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Clear", "[", "n", "]"}], ";"}]}], "Input",
 CellChangeTimes->{{3.5407390737883897`*^9, 3.540739088359223*^9}, {
   3.540739992029157*^9, 3.540740041489399*^9}, {3.5407406597673416`*^9, 
   3.5407406605630484`*^9}, {3.540759923421386*^9, 3.54075992699477*^9}, {
   3.540916666913783*^9, 3.5409167107412057`*^9}, {3.5410306226774387`*^9, 
   3.5410306254230385`*^9}, {3.541270497166806*^9, 3.5412704996322856`*^9}, {
   3.541274869409685*^9, 3.541274944053684*^9}, {3.54127535250709*^9, 
   3.5412753550661683`*^9}, {3.5412756166532516`*^9, 3.541275618503356*^9}, 
   3.541296481310985*^9, {3.5413395075363216`*^9, 3.541339527972351*^9}, {
   3.541854954893719*^9, 3.5418550048764396`*^9}, {3.541901835934777*^9, 
   3.5419018359972024`*^9}, {3.5425556289495316`*^9, 3.542555668948188*^9}, {
   3.5425645375795045`*^9, 3.5425645423687353`*^9}, {3.5425708187683387`*^9, 
   3.542570829282873*^9}, {3.5429947403438854`*^9, 3.5429947719235477`*^9}, {
   3.5429967242128005`*^9, 3.54299673992482*^9}, 3.5429987707942295`*^9, {
   3.54860722625*^9, 3.548607242109375*^9}, {3.5486142270302734`*^9, 
   3.5486142323271484`*^9}, {3.549047229984375*^9, 3.54904723009375*^9}, {
   3.549047281875*^9, 3.5490472825*^9}, {3.5492139240859375`*^9, 
   3.5492139246621094`*^9}, {3.549218988041992*^9, 3.549218988958008*^9}}],

Cell[CellGroupData[{

Cell["Speed of Light", "Subsubtitle"],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"c", "=", "29979245800"}], ";"}], 
  RowBox[{"(*", 
   RowBox[{"cm", "/", "sec"}], "*)"}]}]], "Input"]
}, Open  ]],

Cell[CellGroupData[{

Cell["Planck's Constant", "Subsubtitle"],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"hBar", "=", 
    RowBox[{"1.054", " ", 
     SuperscriptBox["10", 
      RowBox[{"-", "27"}]]}]}], ";"}], 
  RowBox[{"(*", 
   RowBox[{"ergs", " ", "sec"}], "*)"}]}]], "Input"]
}, Open  ]],

Cell[CellGroupData[{

Cell["Boltzmann Constant", "Subsubtitle",
 CellChangeTimes->{{3.456659575640625*^9, 3.456659597484375*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"kB", "=", 
    RowBox[{"1.380658", " ", 
     SuperscriptBox["10", 
      RowBox[{"-", "16"}]]}]}], ";"}], 
  RowBox[{"(*", 
   RowBox[{"erg", "/", "K"}], "*)"}]}]], "Input",
 CellChangeTimes->{{3.4566596001875*^9, 3.456659625953125*^9}, {
  3.45665968075*^9, 3.456659686359375*^9}}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Atomic Parameters", "Subtitle",
 FontColor->RGBColor[0, 0.356863, 0]],

Cell["General Parameters", "Subtitle",
 CellChangeTimes->{{3.45666146840625*^9, 3.45666147296875*^9}}],

Cell["Nuclear Spin", "Subsubtitle",
 CellChangeTimes->{{3.456659575640625*^9, 3.456659597484375*^9}, {
  3.45665970596875*^9, 3.456659711984375*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"i", "=", 
   RowBox[{"N", "[", 
    FractionBox["3", "2"], "]"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.45665971353125*^9, 3.45665971846875*^9}, {
   3.4730885615*^9, 3.473088561953125*^9}, 3.53632801824225*^9, 
   3.5386065281276674`*^9, {3.542513177548191*^9, 3.5425131799194365`*^9}}],

Cell[CellGroupData[{

Cell["Transit Linewidth", "Subsubtitle",
 CellChangeTimes->{{3.45666389359375*^9, 3.456663897328125*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"\[Gamma]T", "=", 
   FractionBox["1", 
    RowBox[{"40.2", " ", 
     SuperscriptBox["10", 
      RowBox[{"-", "9"}]]}]]}], ";"}]], "Input",
 CellChangeTimes->{{3.456663726421875*^9, 3.4566637388125*^9}, {
   3.45666380203125*^9, 3.456663813921875*^9}, 3.4566922269375*^9, 
   3.5366835360235*^9, 3.538523997111776*^9, {3.5492139427109375`*^9, 
   3.5492139432939453`*^9}}],

Cell[CellGroupData[{

Cell[BoxData["\[Gamma]T"], "Input",
 CellChangeTimes->{{3.5370977448424997`*^9, 3.5370977478424997`*^9}}],

Cell[BoxData["2.487562189054726`*^7"], "Output",
 CellChangeTimes->{3.5492150784335938`*^9, 3.549218996876953*^9, 
  3.549221390517578*^9, 3.549223597060547*^9}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Polarization", "Subsubtitle",
 CellChangeTimes->{{3.456663909203125*^9, 3.456663910671875*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"q1", "=", 
   RowBox[{"{", 
    RowBox[{"0", ",", "1", ",", "0"}], "}"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.456663912765625*^9, 3.456663917453125*^9}, 
   3.456663957921875*^9}],

Cell[BoxData[
 RowBox[{
  RowBox[{"q2", "=", 
   RowBox[{"{", 
    RowBox[{"0", ",", "1", ",", "0"}], "}"}]}], ";"}]], "Input",
 CellChangeTimes->{3.4566639629375*^9}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Ground State Parameters", "Subtitle",
 CellChangeTimes->{{3.456659852765625*^9, 3.456659860875*^9}}],

Cell["Ground State S", "Subsubtitle",
 CellChangeTimes->{{3.456659575640625*^9, 3.456659597484375*^9}, {
   3.45665970596875*^9, 3.456659711984375*^9}, {3.45665974584375*^9, 
   3.456659752203125*^9}, 3.45665980753125*^9}],

Cell[BoxData[
 RowBox[{
  RowBox[{"SGnd", "=", 
   RowBox[{"N", "[", 
    FractionBox["1", "2"], "]"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.45665971353125*^9, 3.45665971846875*^9}, {
   3.45665975521875*^9, 3.456659757859375*^9}, 3.4566598095625*^9, {
   3.5425131864715624`*^9, 3.54251319049644*^9}}],

Cell[CellGroupData[{

Cell["Ground State J", "Subsubtitle",
 CellChangeTimes->{{3.456659575640625*^9, 3.456659597484375*^9}, {
  3.45665970596875*^9, 3.456659711984375*^9}, {3.45665974584375*^9, 
  3.456659752203125*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"jGnd", "=", 
   RowBox[{"N", "[", 
    FractionBox["1", "2"], "]"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.45665971353125*^9, 3.45665971846875*^9}, {
  3.45665975521875*^9, 3.456659757859375*^9}, {3.5425131924152765`*^9, 
  3.5425131945057173`*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Ground State F Values", "Subsubtitle",
 CellChangeTimes->{{3.45666290321875*^9, 3.456662908015625*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"fMinGnd", "=", 
   RowBox[{"Abs", "[", 
    RowBox[{"jGnd", "-", "i"}], "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"fMaxGnd", "=", 
   RowBox[{"jGnd", "+", "i"}]}], ";"}]}], "Input",
 CellChangeTimes->{{3.456662861609375*^9, 3.456662890515625*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Ground State Hyperfine A", "Subsubtitle",
 CellChangeTimes->{{3.456659575640625*^9, 3.456659597484375*^9}, {
  3.45665970596875*^9, 3.456659711984375*^9}, {3.456659762*^9, 
  3.45665976646875*^9}, {3.456659909703125*^9, 3.45665991996875*^9}, {
  3.45666013296875*^9, 3.456660152140625*^9}, {3.456660398109375*^9, 
  3.45666040421875*^9}, {3.456660756828125*^9, 3.45666075759375*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"aGnd", "=", 
    RowBox[{"SetPrecision", "[", 
     RowBox[{"3417341305.452145", ",", "20"}], "]"}]}], " ", ";"}], 
  RowBox[{"(*", 
   RowBox[{"From", " ", "Steck", " ", "review"}], "*)"}]}]], "Input",
 CellChangeTimes->{{3.45665971353125*^9, 3.45665971846875*^9}, {
  3.45665976928125*^9, 3.45665977578125*^9}, {3.45665992590625*^9, 
  3.45665992978125*^9}, {3.456659964359375*^9, 3.45665998015625*^9}, {
  3.456660078296875*^9, 3.456660092390625*^9}, {3.456660306609375*^9, 
  3.45666036596875*^9}, {3.456660407359375*^9, 3.45666043703125*^9}, {
  3.4566605581875*^9, 3.4566605719375*^9}, {3.456660680015625*^9, 
  3.456660705609375*^9}, {3.45666076115625*^9, 3.4566607615625*^9}, {
  3.47308857578125*^9, 3.473088584546875*^9}, {3.47309640846875*^9, 
  3.473096424140625*^9}, {3.536335670945375*^9, 3.53633567542975*^9}, {
  3.5399858315742908`*^9, 3.5399858502631235`*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Ground State Hyperfine B", "Subsubtitle",
 CellChangeTimes->{{3.456659575640625*^9, 3.456659597484375*^9}, {
   3.45665970596875*^9, 3.456659711984375*^9}, {3.456659762*^9, 
   3.45665976646875*^9}, {3.456659909703125*^9, 3.45665991996875*^9}, {
   3.45666013296875*^9, 3.456660152140625*^9}, {3.456660398109375*^9, 
   3.45666040421875*^9}, 3.4566605775*^9, {3.45666076390625*^9, 
   3.45666076546875*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"bGnd", "=", "0"}], ";"}]], "Input",
 CellChangeTimes->{{3.45665971353125*^9, 3.45665971846875*^9}, {
  3.45665976928125*^9, 3.45665977578125*^9}, {3.45665992590625*^9, 
  3.45665992978125*^9}, {3.456659964359375*^9, 3.45665998015625*^9}, {
  3.456660078296875*^9, 3.456660092390625*^9}, {3.456660306609375*^9, 
  3.45666036596875*^9}, {3.456660407359375*^9, 3.45666043703125*^9}, {
  3.4566605581875*^9, 3.456660602203125*^9}, {3.45666063840625*^9, 
  3.45666063971875*^9}, {3.45666076803125*^9, 3.456660773828125*^9}}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["\<\
Intermediate State Parameters 5P1/2\
\>", "Subtitle",
 CellChangeTimes->{{3.456659852765625*^9, 3.456659887765625*^9}, {
  3.5431807549355526`*^9, 3.543180758321512*^9}}],

Cell["Intermediate State S", "Subsubtitle",
 CellChangeTimes->{{3.456659575640625*^9, 3.456659597484375*^9}, {
  3.45665970596875*^9, 3.456659711984375*^9}, {3.45665974584375*^9, 
  3.456659752203125*^9}, {3.45665980753125*^9, 3.45665981840625*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"SInt", "=", 
   RowBox[{"N", "[", 
    FractionBox["1", "2"], "]"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.45665971353125*^9, 3.45665971846875*^9}, {
  3.45665975521875*^9, 3.456659757859375*^9}, {3.4566598095625*^9, 
  3.45665982178125*^9}, {3.5425131990610046`*^9, 3.5425132012294464`*^9}}],

Cell[CellGroupData[{

Cell["Intermediate State J", "Subsubtitle",
 CellChangeTimes->{{3.456659575640625*^9, 3.456659597484375*^9}, {
  3.45665970596875*^9, 3.456659711984375*^9}, {3.456659762*^9, 
  3.45665976646875*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"jInt", "=", 
   RowBox[{"N", "[", 
    FractionBox["1", "2"], "]"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.45665971353125*^9, 3.45665971846875*^9}, {
   3.45665976928125*^9, 3.45665977578125*^9}, {3.473088588953125*^9, 
   3.47308858946875*^9}, 3.536333582476625*^9, {3.542513204380707*^9, 
   3.542513206361945*^9}, {3.543180749833208*^9, 3.5431807499424324`*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Intermediate State F Values", "Subsubtitle",
 CellChangeTimes->{{3.45666290321875*^9, 3.456662923203125*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"fMinInt", "=", 
   RowBox[{"Abs", "[", 
    RowBox[{"jInt", "-", "i"}], "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"fMaxInt", "=", 
   RowBox[{"jInt", "+", "i"}]}], ";"}]}], "Input",
 CellChangeTimes->{{3.456662861609375*^9, 3.456662890515625*^9}, {
  3.456662927125*^9, 3.456662936640625*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Intermediate State Linewidth", "Subsubtitle",
 CellChangeTimes->{{3.456659575640625*^9, 3.456659597484375*^9}, {
  3.45665970596875*^9, 3.456659711984375*^9}, {3.456659762*^9, 
  3.45665976646875*^9}, {3.456659909703125*^9, 3.45665991996875*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"\[Gamma]Int", "=", 
    FractionBox["1", 
     RowBox[{"26.2", " ", 
      SuperscriptBox["10", 
       RowBox[{"-", "9"}]]}]]}], ";"}], 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{"Seconds", " ", "From", " ", "Budker"}], ",", "Kimball", ",", " ",
     "DeMille"}], "*)"}]}]], "Input",
 CellChangeTimes->{{3.45665971353125*^9, 3.45665971846875*^9}, {
  3.45665976928125*^9, 3.45665977578125*^9}, {3.45665992590625*^9, 
  3.45665992978125*^9}, {3.456659964359375*^9, 3.45665998015625*^9}, {
  3.456660078296875*^9, 3.456660092390625*^9}, {3.53668310449225*^9, 
  3.536683106007875*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Intermediate State Center of Gravity Energy", "Subsubtitle",
 CellChangeTimes->{{3.456659575640625*^9, 3.456659597484375*^9}, {
  3.45665970596875*^9, 3.456659711984375*^9}, {3.456659762*^9, 
  3.45665976646875*^9}, {3.456659909703125*^9, 3.45665991996875*^9}, {
  3.45666013296875*^9, 3.456660152140625*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"\[Nu]Int", "=", 
   RowBox[{"SetPrecision", "[", 
    RowBox[{"377107463500000", ",", "20"}], "]"}]}], " ", 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"(*", 
   RowBox[{"Hz", " ", "from", " ", "Steck"}], "*)"}]}]}], "Input",
 CellChangeTimes->{{3.45665971353125*^9, 3.45665971846875*^9}, {
   3.45665976928125*^9, 3.45665977578125*^9}, {3.45665992590625*^9, 
   3.45665992978125*^9}, {3.456659964359375*^9, 3.45665998015625*^9}, {
   3.456660078296875*^9, 3.456660092390625*^9}, {3.456660306609375*^9, 
   3.45666036596875*^9}, {3.473088596515625*^9, 3.473088616875*^9}, {
   3.47308883840625*^9, 3.473088843484375*^9}, {3.473096433640625*^9, 
   3.473096439859375*^9}, {3.53633478836725*^9, 3.53633479311725*^9}, 
   3.536607319211*^9, {3.536607447226625*^9, 3.536607498007875*^9}, {
   3.53666646042975*^9, 3.5366664651485*^9}, {3.537893167558856*^9, 
   3.5378931893520937`*^9}, 3.5378943549667425`*^9, 3.5378974743522215`*^9, {
   3.5378985672711487`*^9, 3.53789856825395*^9}, {3.5378998431661897`*^9, 
   3.537899869374236*^9}, 3.5379017571978636`*^9, 3.538606282516619*^9, {
   3.5399857601885653`*^9, 3.5399857692989817`*^9}, {3.543180797673539*^9, 
   3.543180808923662*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Intermediate State Hyperfine A", "Subsubtitle",
 CellChangeTimes->{{3.456659575640625*^9, 3.456659597484375*^9}, {
  3.45665970596875*^9, 3.456659711984375*^9}, {3.456659762*^9, 
  3.45665976646875*^9}, {3.456659909703125*^9, 3.45665991996875*^9}, {
  3.45666013296875*^9, 3.456660152140625*^9}, {3.456660398109375*^9, 
  3.45666040421875*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"aInt", "=", 
    RowBox[{"SetPrecision", "[", 
     RowBox[{"408328000", ",", "20"}], "]"}]}], " ", ";"}], 
  RowBox[{"(*", 
   RowBox[{"Hz", " ", "From", " ", "Steck"}], "*)"}]}]], "Input",
 CellChangeTimes->{{3.45665971353125*^9, 3.45665971846875*^9}, {
  3.45665976928125*^9, 3.45665977578125*^9}, {3.45665992590625*^9, 
  3.45665992978125*^9}, {3.456659964359375*^9, 3.45665998015625*^9}, {
  3.456660078296875*^9, 3.456660092390625*^9}, {3.456660306609375*^9, 
  3.45666036596875*^9}, {3.456660407359375*^9, 3.45666043703125*^9}, {
  3.4566605581875*^9, 3.4566605719375*^9}, {3.473088627078125*^9, 
  3.4730886284375*^9}, {3.473096443796875*^9, 3.473096448328125*^9}, {
  3.5363357313985*^9, 3.536335734726625*^9}, {3.536607079039125*^9, 
  3.53660709955475*^9}, {3.536666472757875*^9, 3.536666476257875*^9}, {
  3.536683146882875*^9, 3.5366831523985*^9}, {3.5372022874831247`*^9, 
  3.5372022919049997`*^9}, {3.539985799703435*^9, 3.5399858030418406`*^9}, {
  3.543180812216001*^9, 3.543180817926882*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Intermediate State Hyperfine B", "Subsubtitle",
 CellChangeTimes->{{3.456659575640625*^9, 3.456659597484375*^9}, {
   3.45665970596875*^9, 3.456659711984375*^9}, {3.456659762*^9, 
   3.45665976646875*^9}, {3.456659909703125*^9, 3.45665991996875*^9}, {
   3.45666013296875*^9, 3.456660152140625*^9}, {3.456660398109375*^9, 
   3.45666040421875*^9}, 3.4566605775*^9}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"bInt", "=", 
    RowBox[{"SetPrecision", "[", 
     RowBox[{"0", ",", "20"}], "]"}]}], " ", ";"}], 
  RowBox[{"(*", 
   RowBox[{"Hz", " ", "From", " ", "Steck"}], "*)"}]}]], "Input",
 CellChangeTimes->{{3.45665971353125*^9, 3.45665971846875*^9}, {
   3.45665976928125*^9, 3.45665977578125*^9}, {3.45665992590625*^9, 
   3.45665992978125*^9}, {3.456659964359375*^9, 3.45665998015625*^9}, {
   3.456660078296875*^9, 3.456660092390625*^9}, {3.456660306609375*^9, 
   3.45666036596875*^9}, {3.456660407359375*^9, 3.45666043703125*^9}, {
   3.4566605581875*^9, 3.456660602203125*^9}, {3.45666063840625*^9, 
   3.45666063971875*^9}, 3.473088631296875*^9, {3.473096452546875*^9, 
   3.473096456171875*^9}, {3.536335742211*^9, 3.536335745961*^9}, {
   3.536607112711*^9, 3.536607116664125*^9}, 3.53666648311725*^9, {
   3.5366831642735*^9, 3.53668317030475*^9}, {3.5372023036862497`*^9, 
   3.5372023085924997`*^9}, {3.539985808704651*^9, 3.5399858160678635`*^9}, {
   3.543180822389483*^9, 3.543180824168282*^9}}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Final State Parameters", "Subtitle",
 CellChangeTimes->{{3.456659852765625*^9, 3.456659887765625*^9}, {
  3.456660101109375*^9, 3.4566601018125*^9}}],

Cell["Final State S", "Subsubtitle",
 CellChangeTimes->{{3.456659575640625*^9, 3.456659597484375*^9}, {
  3.45665970596875*^9, 3.456659711984375*^9}, {3.45665974584375*^9, 
  3.456659752203125*^9}, {3.45665980753125*^9, 3.456659829375*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"SFin", "=", 
   RowBox[{"N", "[", 
    FractionBox["1", "2"], "]"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.45665971353125*^9, 3.45665971846875*^9}, {
   3.45665975521875*^9, 3.456659757859375*^9}, {3.4566598095625*^9, 
   3.45665983434375*^9}, {3.536335431601625*^9, 3.536335431726625*^9}, 
   3.53634248861725*^9, {3.54251321336648*^9, 3.542513215254116*^9}}],

Cell[CellGroupData[{

Cell["Final State J", "Subsubtitle",
 CellChangeTimes->{{3.456659575640625*^9, 3.456659597484375*^9}, {
  3.45665970596875*^9, 3.456659711984375*^9}, {3.456659783546875*^9, 
  3.456659787671875*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"jFin", "=", 
   RowBox[{"N", "[", 
    FractionBox["3", "2"], "]"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.45665971353125*^9, 3.45665971846875*^9}, {
   3.4566597779375*^9, 3.456659780921875*^9}, {3.456661545015625*^9, 
   3.456661545171875*^9}, 3.4569295484212503`*^9, {3.47308863565625*^9, 
   3.473088640640625*^9}, {3.536335439086*^9, 3.536335439226625*^9}, 
   3.5385845502476044`*^9, {3.54251321806217*^9, 3.542513221525437*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Final State F Values", "Subsubtitle",
 CellChangeTimes->{{3.45666290321875*^9, 3.45666294725*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"fMinFin", "=", 
   RowBox[{"Abs", "[", 
    RowBox[{"jFin", "-", "i"}], "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"fMaxFin", "=", 
   RowBox[{"jFin", "+", "i"}]}], ";"}]}], "Input",
 CellChangeTimes->{{3.456662861609375*^9, 3.456662890515625*^9}, {
  3.456662927125*^9, 3.456662957578125*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Final State Linewidth", "Subsubtitle",
 CellChangeTimes->{{3.456659575640625*^9, 3.456659597484375*^9}, {
  3.45665970596875*^9, 3.456659711984375*^9}, {3.456659762*^9, 
  3.45665976646875*^9}, {3.456659909703125*^9, 3.45665991996875*^9}, {
  3.456659988796875*^9, 3.45665998953125*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"\[Gamma]Fin", "=", 
    FractionBox["1", 
     RowBox[{"241.", " ", 
      SuperscriptBox["10", 
       RowBox[{"-", "9"}]]}]]}], ";"}], 
  RowBox[{"(*", 
   RowBox[{"Seconds", " ", "From", " ", "Heavens"}], "*)"}]}]], "Input",
 CellChangeTimes->{{3.45665971353125*^9, 3.45665971846875*^9}, {
  3.45665976928125*^9, 3.45665977578125*^9}, {3.45665992590625*^9, 
  3.45665992978125*^9}, {3.456659964359375*^9, 3.45665998015625*^9}, {
  3.4566600681875*^9, 3.45666007421875*^9}, {3.4566616303125*^9, 
  3.456661631765625*^9}, {3.4569295617493753`*^9, 3.4569295620150003`*^9}, {
  3.45951917796875*^9, 3.45951917865625*^9}, {3.536683236757875*^9, 
  3.536683244726625*^9}, {3.538582219701305*^9, 3.5385822378273416`*^9}}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Energies, Matrix Elements and Transition Probablities", "Subtitle",
 CellChangeTimes->{{3.456681533796875*^9, 3.456681554234375*^9}}],

Cell["Formula for Calculating Energies with HFS", "Subsubtitle",
 CellChangeTimes->{{3.45666156509375*^9, 3.45666157378125*^9}, {
  3.456662570765625*^9, 3.4566625805625*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"\[Kappa]", "[", 
    RowBox[{"f_", ",", "j_"}], "]"}], "=", 
   RowBox[{
    RowBox[{"f", 
     RowBox[{"(", 
      RowBox[{"f", "+", "1"}], ")"}]}], "-", 
    RowBox[{"j", 
     RowBox[{"(", 
      RowBox[{"j", "+", "1"}], ")"}]}], "-", 
    RowBox[{"i", 
     RowBox[{"(", 
      RowBox[{"i", "+", "1"}], ")"}]}]}]}], ";"}]], "Input",
 CellChangeTimes->{{3.45666169846875*^9, 3.456661713*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"\[Delta]E", "[", 
    RowBox[{"a_", ",", "b_", ",", "f_", ",", "j_"}], "]"}], ":=", 
   RowBox[{"If", "[", 
    RowBox[{
     RowBox[{
      RowBox[{"j", "\[NotEqual]", " ", 
       FractionBox["1", "2"]}], "&&", 
      RowBox[{"i", "\[NotEqual]", 
       FractionBox["1", "2"]}]}], ",", 
     RowBox[{
      RowBox[{
       FractionBox["1", "2"], "a", " ", 
       RowBox[{"\[Kappa]", "[", 
        RowBox[{"f", ",", "j"}], "]"}]}], "+", 
      RowBox[{"b", " ", 
       FractionBox[
        RowBox[{
         RowBox[{
          FractionBox["3", "4"], 
          SuperscriptBox[
           RowBox[{"\[Kappa]", "[", 
            RowBox[{"f", ",", "j"}], "]"}], "2"]}], "+", 
         RowBox[{
          FractionBox["3", "4"], 
          RowBox[{"\[Kappa]", "[", 
           RowBox[{"f", ",", "j"}], "]"}]}], "-", 
         RowBox[{"i", 
          RowBox[{"(", 
           RowBox[{"i", "+", "1"}], ")"}], "j", 
          RowBox[{"(", 
           RowBox[{"j", "+", "1"}], ")"}]}]}], 
        RowBox[{"2", " ", "i", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"2", "i"}], "-", "1"}], ")"}], "j", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"2", "j"}], "-", "1"}], ")"}]}]]}]}], ",", 
     RowBox[{
      FractionBox["1", "2"], "a", " ", 
      RowBox[{"\[Kappa]", "[", 
       RowBox[{"f", ",", "j"}], "]"}]}]}], "]"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.456661820734375*^9, 3.45666183171875*^9}, {
   3.456662393875*^9, 3.456662397109375*^9}, {3.456683422046875*^9, 
   3.456683431546875*^9}, 3.460223934290625*^9}],

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"\[Nu]GndToInt", "[", 
    RowBox[{"f_", ",", "fPrime_"}], "]"}], "=", 
   RowBox[{"\[Nu]Int", "-", 
    RowBox[{"\[Delta]E", "[", 
     RowBox[{"aGnd", ",", "bGnd", ",", "f", ",", "jGnd"}], "]"}], "+", 
    RowBox[{"\[Delta]E", "[", 
     RowBox[{"aInt", ",", "bInt", ",", "fPrime", ",", "jInt"}], "]"}]}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"\[Nu]GndToFin", "[", 
   RowBox[{"f_", ",", "fDPrime_"}], "]"}], ":=", 
  RowBox[{"2", " ", 
   RowBox[{"Switch", "[", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{"f", ",", "fDPrime"}], "}"}], ",", "\[IndentingNewLine]", 
     RowBox[{"{", 
      RowBox[{"2.", ",", "0."}], "}"}], ",", "385240094977700", ",", 
     "\[IndentingNewLine]", 
     RowBox[{"{", 
      RowBox[{"2.", ",", "1."}], "}"}], ",", "385240101727600", ",", 
     "\[IndentingNewLine]", 
     RowBox[{"{", 
      RowBox[{"2.", ",", "2."}], "}"}], ",", "385240115692600", ",", 
     "\[IndentingNewLine]", 
     RowBox[{"{", 
      RowBox[{"2.", ",", "3."}], "}"}], ",", "385240137803800", ",", 
     "\[IndentingNewLine]", 
     RowBox[{"{", 
      RowBox[{"1.", ",", "0."}], "}"}], ",", "385243512320300", ",", 
     "\[IndentingNewLine]", 
     RowBox[{"{", 
      RowBox[{"1.", ",", "1."}], "}"}], ",", "385243519070200", ",", 
     "\[IndentingNewLine]", 
     RowBox[{"{", 
      RowBox[{"1.", ",", "2."}], "}"}], ",", "385243533031800", ",", 
     "\[IndentingNewLine]", 
     RowBox[{"{", 
      RowBox[{"1.", ",", "3."}], "}"}], ",", "385243555144500"}], 
    "]"}]}]}]}], "Input",
 CellChangeTimes->{{3.45686233334375*^9, 3.4568623374375*^9}, 
   3.4572661415462503`*^9, {3.47308885534375*^9, 3.473088860125*^9}, {
   3.4730964853125*^9, 3.473096517875*^9}, {3.473097064765625*^9, 
   3.47309707559375*^9}, {3.536683396664125*^9, 3.53668342911725*^9}, {
   3.5385822740464134`*^9, 3.5385822762455387`*^9}, {3.538582433139513*^9, 
   3.538582435955674*^9}, {3.53858286930646*^9, 3.538583146448312*^9}, {
   3.5386057393735533`*^9, 3.5386058876180325`*^9}, {3.5386065596574707`*^9, 
   3.538606679034299*^9}, {3.538688905996657*^9, 3.5386889171506767`*^9}, {
   3.542513225909121*^9, 3.5425132352068996`*^9}}],

Cell[CellGroupData[{

Cell["Reduced Matrix Elements in J Basis", "Subsubtitle"],

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"dRedJGndToInt", "=", "1"}], ";"}], 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{
     RowBox[{"(", "4.51", ")"}], "1.28", " ", 
     SuperscriptBox["10", "6"]}], " ", 
    RowBox[{"(*", 
     RowBox[{"Hz", "/", 
      RowBox[{"(", 
       RowBox[{"V", "/", "cm"}], ")"}]}], "*)"}], ";"}], 
   "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"dRedJIntToFin", "=", "1"}], ";"}], 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{
     FractionBox["1", 
      RowBox[{"e", " ", "a0"}]], 
     RowBox[{"\[Sqrt]", 
      RowBox[{"(", 
       RowBox[{
        FractionBox[
         RowBox[{
          RowBox[{"(", 
           RowBox[{
            RowBox[{"2", " ", 
             FractionBox["1", "2"]}], "+", "1"}], " ", ")"}], "3", " ", 
          "hBar", " ", 
          SuperscriptBox["c", "3"]}], 
         RowBox[{"4", " ", 
          SuperscriptBox[
           RowBox[{"(", 
            RowBox[{"2", " ", "\[Pi]", " ", "\[Nu]6P12To8S"}], ")"}], "3"]}]],
         "2.04", " ", 
        SuperscriptBox["10", "6"]}], ")"}]}], "1.28", " ", 
     SuperscriptBox["10", "6"]}], " ", 
    RowBox[{"(*", 
     RowBox[{"Hz", "/", 
      RowBox[{"(", 
       RowBox[{"V", "/", "cm"}], ")"}]}], "*)"}], ";"}], "*)"}]}]}], "Input",
 CellChangeTimes->{{3.4566056364375*^9, 3.4566056744375*^9}, {
  3.456662686453125*^9, 3.456662724453125*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Reduced Matrix Elements in F Basis", "Subsubtitle",
 CellChangeTimes->{3.4566632656875*^9}],

Cell[BoxData[
 RowBox[{"Do", "[", "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{"Do", "[", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{
      RowBox[{
       SubscriptBox["dRedGndToInt", 
        RowBox[{"f", ",", "fPrime"}]], "=", 
       RowBox[{"\[Sqrt]", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{
            RowBox[{"2", " ", "fPrime"}], "+", "1"}], ")"}], 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"2", "f"}], "+", "1"}], ")"}], 
          SuperscriptBox[
           RowBox[{"SixJSymbol", "[", 
            RowBox[{
             RowBox[{"{", 
              RowBox[{"jInt", ",", "fPrime", ",", "i"}], "}"}], ",", 
             RowBox[{"{", 
              RowBox[{"f", ",", "jGnd", ",", "1"}], "}"}]}], "]"}], "2"], 
          SuperscriptBox["dRedJGndToInt", "2"]}], ")"}]}]}], ";"}], ",", 
     RowBox[{"{", 
      RowBox[{"fPrime", ",", "fMinInt", ",", "fMaxInt"}], "}"}]}], "]"}], ",",
    "\[IndentingNewLine]", 
   RowBox[{"{", 
    RowBox[{"f", ",", "fMinGnd", ",", "fMaxGnd"}], "}"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.456662795671875*^9, 3.456662835546875*^9}, {
  3.456663002453125*^9, 3.456663020171875*^9}, {3.5372063527331247`*^9, 
  3.5372063614518747`*^9}}],

Cell[BoxData[
 RowBox[{"Do", "[", "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{"Do", "[", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{
      SubscriptBox["dRedIntToFin", 
       RowBox[{"fPrime", ",", "fDPrime"}]], "=", 
      RowBox[{"\[Sqrt]", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{
           RowBox[{"2", " ", "fDPrime"}], "+", "1"}], ")"}], 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"2", "fPrime"}], "+", "1"}], ")"}], 
         SuperscriptBox[
          RowBox[{"SixJSymbol", "[", 
           RowBox[{
            RowBox[{"{", 
             RowBox[{"jFin", ",", "fDPrime", ",", "i"}], "}"}], ",", 
            RowBox[{"{", 
             RowBox[{"fPrime", ",", "jInt", ",", "1"}], "}"}]}], "]"}], "2"], 
         SuperscriptBox["dRedJIntToFin", "2"]}], ")"}]}]}], ",", 
     RowBox[{"{", 
      RowBox[{"fDPrime", ",", "fMinFin", ",", "fMaxFin"}], "}"}]}], "]"}], 
   ",", 
   RowBox[{"{", 
    RowBox[{"fPrime", ",", "fMinInt", ",", "fMaxInt"}], "}"}]}], 
  "]"}]], "Input",
 CellChangeTimes->{{3.456663051265625*^9, 3.456663093375*^9}, {
   3.456663125046875*^9, 3.456663154296875*^9}, 3.456663192609375*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Matrix Elements in F Basis", "Subsubtitle",
 CellChangeTimes->{{3.456663276328125*^9, 3.456663278140625*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"Do", "[", "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Do", "[", "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{"Do", "[", "\[IndentingNewLine]", 
       RowBox[{
        RowBox[{"Do", "[", "\[IndentingNewLine]", 
         RowBox[{
          RowBox[{
           SubscriptBox["dGndToInt", 
            RowBox[{"f", ",", "fPrime", ",", "m", ",", "mPrime"}]], "=", 
           RowBox[{
            SuperscriptBox[
             RowBox[{"(", 
              RowBox[{"-", "1"}], ")"}], 
             RowBox[{"fPrime", "-", "mPrime"}]], 
            RowBox[{
             UnderoverscriptBox["\[Sum]", 
              RowBox[{"n", "=", "1"}], "3"], 
             RowBox[{
              RowBox[{"(", 
               RowBox[{
                RowBox[{
                "q1", "\[LeftDoubleBracket]", "n", "\[RightDoubleBracket]"}], 
                RowBox[{"ThreeJSymbol", "[", 
                 RowBox[{
                  RowBox[{"{", 
                   RowBox[{"fPrime", ",", 
                    RowBox[{"-", "mPrime"}]}], "}"}], ",", 
                  RowBox[{"{", 
                   RowBox[{"1", ",", 
                    RowBox[{"n", "-", "2"}]}], "}"}], ",", 
                  RowBox[{"{", 
                   RowBox[{"f", ",", "m"}], "}"}]}], "]"}]}], ")"}], 
              SubscriptBox["dRedGndToInt", 
               RowBox[{"f", ",", "fPrime"}]]}]}]}]}], ",", 
          "\[IndentingNewLine]", 
          RowBox[{"{", 
           RowBox[{"m", ",", 
            RowBox[{"-", "f"}], ",", "f"}], "}"}]}], "]"}], ",", 
        "\[IndentingNewLine]", 
        RowBox[{"{", 
         RowBox[{"mPrime", ",", 
          RowBox[{"-", "fPrime"}], ",", "fPrime"}], "}"}]}], "]"}], ",", 
      "\[IndentingNewLine]", 
      RowBox[{"{", 
       RowBox[{"fPrime", ",", "fMinInt", ",", "fMaxInt"}], "}"}]}], "]"}], 
    ",", "\[IndentingNewLine]", 
    RowBox[{"{", 
     RowBox[{"f", ",", "fMinGnd", ",", "fMaxGnd"}], "}"}]}], "]"}], 
  ";"}]], "Input",
 CellChangeTimes->{{3.45666329196875*^9, 3.4566633891875*^9}, {3.456663993*^9,
   3.45666408128125*^9}, {3.45666411909375*^9, 3.456664121515625*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"Do", "[", "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Do", "[", "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{"Do", "[", "\[IndentingNewLine]", 
       RowBox[{
        RowBox[{"Do", "[", "\[IndentingNewLine]", 
         RowBox[{
          RowBox[{
           SubscriptBox["dIntToFin", 
            RowBox[{
            "fPrime", ",", "fDPrime", ",", "mPrime", ",", "mDPrime"}]], "=", 
           RowBox[{
            SuperscriptBox[
             RowBox[{"(", 
              RowBox[{"-", "1"}], ")"}], 
             RowBox[{"fDPrime", "-", "mDPrime"}]], 
            RowBox[{
             UnderoverscriptBox["\[Sum]", 
              RowBox[{"n", "=", "1"}], "3"], 
             RowBox[{
              RowBox[{"(", 
               RowBox[{
                RowBox[{
                "q2", "\[LeftDoubleBracket]", "n", "\[RightDoubleBracket]"}], 
                RowBox[{"ThreeJSymbol", "[", 
                 RowBox[{
                  RowBox[{"{", 
                   RowBox[{"fDPrime", ",", 
                    RowBox[{"-", "mDPrime"}]}], "}"}], ",", 
                  RowBox[{"{", 
                   RowBox[{"1", ",", 
                    RowBox[{"n", "-", "2"}]}], "}"}], ",", 
                  RowBox[{"{", 
                   RowBox[{"fPrime", ",", "mPrime"}], "}"}]}], "]"}]}], ")"}], 
              SubscriptBox["dRedIntToFin", 
               RowBox[{"fPrime", ",", "fDPrime"}]]}]}]}]}], ",", 
          "\[IndentingNewLine]", 
          RowBox[{"{", 
           RowBox[{"mPrime", ",", 
            RowBox[{"-", "fPrime"}], ",", "fPrime"}], "}"}]}], "]"}], ",", 
        "\[IndentingNewLine]", 
        RowBox[{"{", 
         RowBox[{"mDPrime", ",", 
          RowBox[{"-", "fDPrime"}], ",", "fDPrime"}], "}"}]}], "]"}], ",", 
      "\[IndentingNewLine]", 
      RowBox[{"{", 
       RowBox[{"fDPrime", ",", "fMinFin", ",", "fMaxFin"}], "}"}]}], "]"}], 
    ",", "\[IndentingNewLine]", 
    RowBox[{"{", 
     RowBox[{"fPrime", ",", "fMinInt", ",", "fMaxInt"}], "}"}]}], "]"}], 
  ";"}]], "Input",
 CellChangeTimes->{{3.45666329196875*^9, 3.456663489203125*^9}, {
  3.456664105703125*^9, 3.45666415221875*^9}}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Doppler Distribution Parameters", "Subtitle",
 FontColor->RGBColor[0, 0.356863, 0]],

Cell[BoxData[{
 RowBox[{
  RowBox[{"amu", "=", 
   RowBox[{
    RowBox[{"(", 
     RowBox[{"1.66053886", " ", 
      SuperscriptBox["10", 
       RowBox[{"-", "27"}]]}], ")"}], " ", 
    SuperscriptBox["10", "3"]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"mRb", "=", 
   RowBox[{"86.909180520", " ", "amu"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"temp", "=", "323"}], ";"}]}], "Input",
 CellChangeTimes->{{3.456606503203125*^9, 3.4566065051875*^9}, {
   3.456664763890625*^9, 3.456664790515625*^9}, {3.473097352640625*^9, 
   3.473097368078125*^9}, {3.536605457211*^9, 3.53660545749225*^9}, 
   3.5366060920235*^9, {3.536606131101625*^9, 3.536606151414125*^9}, 
   3.53666517880475*^9, {3.5402470216155534`*^9, 3.540247090642502*^9}, {
   3.540668358386324*^9, 3.5406683636016226`*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"\[Alpha]", "=", 
  RowBox[{"\[Sqrt]", 
   RowBox[{"(", 
    FractionBox[
     RowBox[{"2", " ", "kB", " ", "temp"}], "mRb"], ")"}]}]}]], "Input",
 CellChangeTimes->{{3.45666479665625*^9, 3.456664796984375*^9}, {
  3.536606154070375*^9, 3.536606154336*^9}}],

Cell[BoxData["24860.048218994714`"], "Output",
 CellChangeTimes->{3.5492150801689453`*^9, 3.5492189982890625`*^9, 
  3.5492213918789062`*^9, 3.5492235983359375`*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{"dopplerDist", "[", "v_", "]"}], ":=", 
  RowBox[{"Exp", "[", 
   RowBox[{"-", 
    FractionBox[
     SuperscriptBox["v", "2"], 
     SuperscriptBox["\[Alpha]", "2"]]}], "]"}]}]], "Input"]
}, Open  ]],

Cell[CellGroupData[{

Cell["Comb Parameters", "Subtitle",
 FontColor->RGBColor[0, 0.356863, 0]],

Cell["Offset Frequency", "Subsubtitle"],

Cell[BoxData[{
 RowBox[{
  RowBox[{"f0", "=", 
   RowBox[{
    RowBox[{"-", "40"}], "*", 
    SuperscriptBox["10", "6"]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"scanRange", "=", "90"}], ";"}]}], "Input",
 CellChangeTimes->{{3.536606229336*^9, 3.536606229476625*^9}, 
   3.536609033945375*^9, {3.53667308924225*^9, 3.536673089507875*^9}, {
   3.536679694632875*^9, 3.5366796956485*^9}, {3.536683666961*^9, 
   3.5366836672735*^9}, 3.5370104677487497`*^9, 3.5370325434049997`*^9, {
   3.539116492685813*^9, 3.5391164930290136`*^9}, 3.539127919839781*^9, {
   3.539999557079931*^9, 3.539999600291931*^9}, {3.540042597806141*^9, 
   3.5400425992725315`*^9}, 3.5400659159530563`*^9, {3.5425132689347486`*^9, 
   3.542513278107725*^9}, {3.543180774112254*^9, 3.5431807871879873`*^9}}],

Cell[CellGroupData[{

Cell["Rep Rate Range", "Subsubtitle",
 CellChangeTimes->{{3.45666493915625*^9, 3.456664941796875*^9}}],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"fRepMin", "=", "250000000.0"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{"fRepMax", "=", 
  RowBox[{"fRepMin", "+", "scanRange"}]}]}], "Input",
 CellChangeTimes->{{3.456664947171875*^9, 3.4566649499375*^9}, {
   3.456665042671875*^9, 3.45666507275*^9}, {3.456665103265625*^9, 
   3.45666510571875*^9}, {3.456665200171875*^9, 3.45666520296875*^9}, {
   3.456691746078125*^9, 3.45669175240625*^9}, {3.536342589414125*^9, 
   3.536342618445375*^9}, {3.5366062423985*^9, 3.536606269101625*^9}, {
   3.536673100086*^9, 3.536673151445375*^9}, {3.536679707945375*^9, 
   3.536679725882875*^9}, {3.536683688101625*^9, 3.53668372074225*^9}, {
   3.5369246547956247`*^9, 3.5369246550456247`*^9}, {3.5372003428581247`*^9, 
   3.5372003561081247`*^9}, {3.5377186875620003`*^9, 
   3.5377186920307503`*^9}, {3.5377228732026253`*^9, 
   3.5377228764838753`*^9}, {3.5377945896870003`*^9, 
   3.5377945969526253`*^9}, {3.5378884679682503`*^9, 
   3.5378884728745003`*^9}, {3.5378886446870003`*^9, 
   3.5378886463745003`*^9}, {3.5378933335587473`*^9, 
   3.5378933356803513`*^9}, {3.537893370718013*^9, 3.5378933755228214`*^9}, {
   3.5379017917519245`*^9, 3.537901792812726*^9}, {3.5385833273276577`*^9, 
   3.538583332126932*^9}, {3.539116509315442*^9, 3.539116540281497*^9}, {
   3.539999589730731*^9, 3.539999598310731*^9}, 3.539999698417289*^9, 
   3.540065826867961*^9, 3.540065927064692*^9, 3.540077153730686*^9, {
   3.5404089194488587`*^9, 3.540408922466031*^9}}],

Cell[BoxData["2.5000009`*^8"], "Output",
 CellChangeTimes->{3.5492150802626953`*^9, 3.5492189983828125`*^9, 
  3.5492213919472656`*^9, 3.5492235984140625`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"fRepAvg", "=", 
  FractionBox[
   RowBox[{"fRepMin", "+", "fRepMax"}], "2"]}]], "Input",
 CellChangeTimes->{{3.456665192578125*^9, 3.456665212796875*^9}, 
   3.5400659248425646`*^9}],

Cell[BoxData["2.50000045`*^8"], "Output",
 CellChangeTimes->{3.5492150803554688`*^9, 3.5492189984492188`*^9, 
  3.549221391998047*^9, 3.5492235984697266`*^9}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Comb Modes", "Subsubtitle",
 CellChangeTimes->{{3.4566813895*^9, 3.45668139175*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"n1NomA", "=", 
   RowBox[{"Round", "[", 
    FractionBox[
     RowBox[{"\[Nu]GndToInt", "[", 
      RowBox[{"fMinGnd", ",", "fMinInt"}], "]"}], "fRepAvg"], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"n1NomB", "=", 
   RowBox[{"Round", "[", 
    FractionBox[
     RowBox[{"\[Nu]GndToInt", "[", 
      RowBox[{"fMaxGnd", ",", "fMaxInt"}], "]"}], "fRepAvg"], "]"}]}], 
  ";"}]}], "Input",
 CellChangeTimes->{{3.456665147765625*^9, 3.456665238296875*^9}, {
  3.5404090137802544`*^9, 3.540409014676306*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"n2NomA", "=", 
   RowBox[{"Round", "[", 
    RowBox[{
     FractionBox["1", "fRepAvg"], 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"\[Nu]GndToFin", "[", 
        RowBox[{"fMinGnd", ",", "fMinFin"}], "]"}], "-", 
       RowBox[{"\[Nu]GndToInt", "[", 
        RowBox[{"fMinGnd", ",", "fMinInt"}], "]"}]}], ")"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"n2NomB", "=", 
   RowBox[{"Round", "[", 
    RowBox[{
     FractionBox["1", "fRepAvg"], 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"\[Nu]GndToFin", "[", 
        RowBox[{"fMaxGnd", ",", "fMaxFin"}], "]"}], "-", 
       RowBox[{"\[Nu]GndToInt", "[", 
        RowBox[{"fMaxGnd", ",", "fMaxInt"}], "]"}]}], ")"}]}], "]"}]}], 
  ";"}]}], "Input",
 CellChangeTimes->{{3.456665147765625*^9, 3.4566653128125*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"n1List", "=", 
  RowBox[{"Union", "[", 
   RowBox[{"Join", "[", 
    RowBox[{
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"n1NomA", "+", "\[Delta]n"}], ",", 
       RowBox[{"{", 
        RowBox[{"\[Delta]n", ",", 
         RowBox[{"-", "3"}], ",", "3"}], "}"}]}], "]"}], ",", 
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"n1NomB", "+", "\[Delta]n"}], ",", 
       RowBox[{"{", 
        RowBox[{"\[Delta]n", ",", 
         RowBox[{"-", "3"}], ",", "3"}], "}"}]}], "]"}]}], "]"}], 
   "]"}]}]], "Input",
 CellChangeTimes->{{3.456665383359375*^9, 3.45666546196875*^9}, {
  3.456665618125*^9, 3.456665620609375*^9}, {3.5372066701081247`*^9, 
  3.5372066730768747`*^9}, {3.5377187815151253`*^9, 3.5377187956088753`*^9}, {
  3.5377946040932503`*^9, 3.5377946077338753`*^9}, {3.5377949385776253`*^9, 
  3.5377949450151253`*^9}, {3.537893207650926*^9, 3.537893217432143*^9}, {
  3.5385234618671618`*^9, 3.53852346744048*^9}, {3.538583211301021*^9, 
  3.5385832213905983`*^9}, {3.5403169380378237`*^9, 3.540316942063656*^9}, {
  3.5404089486965313`*^9, 3.540408953899829*^9}, {3.540408989700877*^9, 
  3.540408998262367*^9}, {3.540442631448593*^9, 3.540442639810139*^9}, {
  3.5425131013095245`*^9, 3.5425131054592047`*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
  "1508418", ",", "1508419", ",", "1508420", ",", "1508421", ",", "1508422", 
   ",", "1508423", ",", "1508424", ",", "1508442", ",", "1508443", ",", 
   "1508444", ",", "1508445", ",", "1508446", ",", "1508447", ",", 
   "1508448"}], "}"}]], "Output",
 CellChangeTimes->{3.5492150804277344`*^9, 3.549218998551758*^9, 
  3.5492213920722656`*^9, 3.549223598544922*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"n2List", "=", 
  RowBox[{"Union", "[", 
   RowBox[{"Join", "[", 
    RowBox[{
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"n2NomA", "+", "\[Delta]n"}], ",", 
       RowBox[{"{", 
        RowBox[{"\[Delta]n", ",", 
         RowBox[{"-", "3"}], ",", "3"}], "}"}]}], "]"}], ",", 
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"n2NomB", "+", "\[Delta]n"}], ",", 
       RowBox[{"{", 
        RowBox[{"\[Delta]n", ",", 
         RowBox[{"-", "3"}], ",", "3"}], "}"}]}], "]"}]}], "]"}], 
   "]"}]}]], "Input",
 CellChangeTimes->{{3.456665383359375*^9, 3.45666546196875*^9}, {
  3.456665585125*^9, 3.45666558903125*^9}, {3.456665624328125*^9, 
  3.456665626734375*^9}, {3.5372066753268747`*^9, 3.5372066804518747`*^9}, {
  3.5377187866401253`*^9, 3.5377187978745003`*^9}, {3.5377946088901253`*^9, 
  3.5377946116245003`*^9}, {3.5377949469370003`*^9, 3.5377949552338753`*^9}, {
  3.5378932198813477`*^9, 3.537893225029357*^9}, {3.538523470201638*^9, 
  3.5385234737538414`*^9}, {3.538583224257762*^9, 3.5385832272179317`*^9}, {
  3.54031694356164*^9, 3.540316948430088*^9}, {3.5404089555279226`*^9, 
  3.5404089605492096`*^9}, {3.5404089910619545`*^9, 3.540408996635274*^9}, {
  3.5404426567672305`*^9, 3.5404426647075796`*^9}, {3.5425131072220383`*^9, 
  3.542513113384157*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
  "1573497", ",", "1573498", ",", "1573499", ",", "1573500", ",", "1573501", 
   ",", "1573502", ",", "1573503", ",", "1573504", ",", "1573505", ",", 
   "1573506"}], "}"}]], "Output",
 CellChangeTimes->{3.5492150804804688`*^9, 3.5492189986289062`*^9, 
  3.549221392125*^9, 3.549223598600586*^9}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Comb Equation", "Subsubtitle",
 CellChangeTimes->{{3.456681410328125*^9, 3.456681413015625*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"\[Omega]L", "[", 
    RowBox[{"n_", ",", "fR_", ",", "f0_"}], "]"}], ":=", 
   RowBox[{"2", " ", "\[Pi]", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"n", " ", "fR"}], "+", "f0"}], ")"}]}]}], ";"}]], "Input",
 CellChangeTimes->{{3.456678515171875*^9, 3.45667856759375*^9}, {
   3.4566795630625*^9, 3.45667956609375*^9}, 3.54904762865625*^9}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Integrate 2-photon Lineshape over Doppler Distribution", "Subtitle",
 FontColor->RGBColor[0, 0.356863, 0]],

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"transProb", "[", 
     RowBox[{"\[Omega]1_", ",", "\[Omega]2_", ",", "v_"}], "]"}], "=", 
    RowBox[{
     UnderoverscriptBox["\[Sum]", 
      RowBox[{"fDPrime", "=", "fMinFin"}], "fMaxFin"], 
     RowBox[{
      UnderoverscriptBox["\[Sum]", 
       RowBox[{"f", "=", "fMinGnd"}], "fMaxGnd"], 
      RowBox[{"(", 
       RowBox[{
        RowBox[{
         FractionBox[
          SuperscriptBox[
           RowBox[{"(", 
            RowBox[{"\[Gamma]Fin", "+", "\[Gamma]T"}], ")"}], "2"], "4"], "/", 
         RowBox[{"(", 
          RowBox[{
           SuperscriptBox[
            RowBox[{"(", 
             RowBox[{
              RowBox[{"2", " ", "\[Pi]", " ", 
               RowBox[{"\[Nu]GndToFin", "[", 
                RowBox[{"f", ",", "fDPrime"}], "]"}]}], "-", 
              RowBox[{
               RowBox[{"(", "\[Omega]1", ")"}], 
               RowBox[{"(", 
                RowBox[{"1", "+", 
                 FractionBox["v", "c"]}], ")"}]}], "-", 
              RowBox[{
               RowBox[{"(", "\[Omega]2", ")"}], 
               RowBox[{"(", 
                RowBox[{"1", "-", 
                 FractionBox["v", "c"]}], ")"}]}]}], ")"}], "2"], "+", 
           FractionBox[
            SuperscriptBox[
             RowBox[{"(", 
              RowBox[{"\[Gamma]Fin", "+", "\[Gamma]T"}], ")"}], "2"], "4"]}], 
          ")"}]}], "\[IndentingNewLine]", 
        FractionBox["1", 
         RowBox[{
          RowBox[{"2", "f"}], "+", "1"}]], 
        RowBox[{
         UnderoverscriptBox["\[Sum]", 
          RowBox[{"m", "=", 
           RowBox[{"-", "f"}]}], "f"], 
         RowBox[{
          UnderoverscriptBox["\[Sum]", 
           RowBox[{"mDPrime", "=", 
            RowBox[{"-", "fDPrime"}]}], "fDPrime"], 
          SuperscriptBox[
           RowBox[{"Abs", "[", 
            RowBox[{
             UnderoverscriptBox["\[Sum]", 
              RowBox[{"fPrime", "=", "fMinInt"}], "fMaxInt"], 
             RowBox[{
              UnderoverscriptBox["\[Sum]", 
               RowBox[{"mPrime", "=", 
                RowBox[{"-", "fPrime"}]}], "fPrime"], "\[IndentingNewLine]", 
              RowBox[{
               RowBox[{"(", 
                RowBox[{
                 FractionBox[
                  RowBox[{"\[Gamma]Int", "+", "\[Gamma]T"}], "2"], "/", 
                 RowBox[{"(", 
                  RowBox[{
                   FractionBox[
                    RowBox[{"\[ImaginaryI]", " ", 
                    RowBox[{"(", 
                    RowBox[{"\[Gamma]Int", "+", "\[Gamma]T"}], ")"}]}], "2"], 
                   "+", 
                   RowBox[{"(", 
                    RowBox[{
                    RowBox[{
                    RowBox[{"(", "\[Omega]1", ")"}], 
                    RowBox[{"(", 
                    RowBox[{"1", "+", 
                    FractionBox["v", "c"]}], ")"}]}], "-", " ", 
                    RowBox[{"2", " ", "\[Pi]", " ", 
                    RowBox[{"\[Nu]GndToInt", "[", 
                    RowBox[{"f", ",", "fPrime"}], "]"}]}]}], ")"}]}], ")"}]}],
                 ")"}], 
               SubscriptBox["dGndToInt", 
                RowBox[{"f", ",", "fPrime", ",", "m", ",", "mPrime"}]], 
               SubscriptBox["dIntToFin", 
                RowBox[{
                "fPrime", ",", "fDPrime", ",", "mPrime", ",", 
                 "mDPrime"}]]}]}]}], "]"}], "2"]}]}]}], ")"}]}]}]}], ";"}], 
  "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"fnToIntegrate", "[", 
     RowBox[{"fR_", ",", "fOffset_", ",", "v_"}], "]"}], "=", 
    RowBox[{
     UnderoverscriptBox["\[Sum]", 
      RowBox[{"n1", "=", "1"}], 
      RowBox[{"Length", "[", "n1List", "]"}]], 
     RowBox[{
      UnderoverscriptBox["\[Sum]", 
       RowBox[{"n2", "=", "1"}], 
       RowBox[{"Length", "[", "n2List", "]"}]], 
      RowBox[{
       RowBox[{"dopplerDist", "[", "v", "]"}], 
       RowBox[{"transProb", "[", 
        RowBox[{
         RowBox[{"\[Omega]L", "[", 
          RowBox[{
           RowBox[{
           "n1List", "\[LeftDoubleBracket]", "n1", "\[RightDoubleBracket]"}], 
           ",", "fR", ",", "fOffset"}], "]"}], ",", 
         RowBox[{"\[Omega]L", "[", 
          RowBox[{
           RowBox[{
           "n2List", "\[LeftDoubleBracket]", "n2", "\[RightDoubleBracket]"}], 
           ",", "fR", ",", "fOffset"}], "]"}], ",", "v"}], "]"}]}]}]}]}], 
   ";"}], "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"vIntSpec", "[", 
    RowBox[{
     RowBox[{"fRep_", "?", "NumericQ"}], ",", 
     RowBox[{"fO_", "?", "NumericQ"}]}], "]"}], ":=", 
   RowBox[{"NIntegrate", "[", 
    RowBox[{
     RowBox[{"fnToIntegrate", "[", 
      RowBox[{"fRep", ",", "fO", ",", "v"}], "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"v", ",", 
       RowBox[{"-", "60000"}], ",", "60000"}], "}"}], ",", 
     RowBox[{"MinRecursion", "\[Rule]", "3"}], ",", 
     RowBox[{"MaxRecursion", "\[Rule]", "10"}], ",", 
     RowBox[{"AccuracyGoal", "\[Rule]", "4"}], ",", 
     RowBox[{"PrecisionGoal", "\[Rule]", "2"}]}], "]"}]}], ";"}]}], "Input",
 CellChangeTimes->{{3.538523548469115*^9, 3.5385235524603434`*^9}, {
  3.5385832601238136`*^9, 3.538583261230877*^9}, {3.54077684144771*^9, 
  3.5407768433509836`*^9}, {3.5407781956265917`*^9, 3.5407781988247147`*^9}, {
  3.540778342428238*^9, 3.540778354378298*^9}, {3.54090952236373*^9, 
  3.5409095316782627`*^9}, {3.5409167735568705`*^9, 3.5409167781128006`*^9}, {
  3.5409170292350388`*^9, 3.5409170310449286`*^9}, {3.5409218751219053`*^9, 
  3.5409218761359835`*^9}, {3.54127538747349*^9, 3.5412754413861923`*^9}, {
  3.5412754775292873`*^9, 3.5412754971543236`*^9}, {3.541854908873424*^9, 
  3.541854930573163*^9}, {3.541901821732953*^9, 3.5419018296297913`*^9}, {
  3.5486142020615234`*^9, 3.5486142102021484`*^9}, {3.5486152438740234`*^9, 
  3.5486152532802734`*^9}, {3.5486153189365234`*^9, 3.5486153292021484`*^9}, {
  3.5486156190771484`*^9, 3.5486156448271484`*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Calculation", "Subtitle",
 CellChangeTimes->{{3.456681446140625*^9, 3.45668144753125*^9}, {
  3.5379007941050625`*^9, 3.537900798020669*^9}},
 FontColor->RGBColor[0, 0.356863, 0]],

Cell[BoxData[
 RowBox[{"(*", 
  RowBox[{
   RowBox[{"\[Delta]f", "=", "0.25"}], ";", " ", 
   RowBox[{"(*", 
    RowBox[{"step", " ", "size", " ", "in", " ", "rep", " ", "rate"}], "*)"}],
    "\[IndentingNewLine]", 
   RowBox[{"fr", "=", "250000000"}], ";", "\[IndentingNewLine]", 
   RowBox[{"fRepList", "=", 
    RowBox[{"{", "}"}]}], ";", "\[IndentingNewLine]", 
   RowBox[{"specList", "=", 
    RowBox[{"{", "}"}]}], ";", "\[IndentingNewLine]", 
   RowBox[{"(*", 
    RowBox[{"Form", " ", "the", " ", "fRepList"}], "*)"}], 
   "\[IndentingNewLine]", 
   RowBox[{"While", " ", "[", 
    RowBox[{
     RowBox[{"fr", "<", "250000090"}], ",", "\[IndentingNewLine]", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"fRepList", "=", 
        RowBox[{"Append", "[", 
         RowBox[{"fRepList", ",", "fr"}], "]"}]}], ";", "\[IndentingNewLine]", 
       RowBox[{"fr", "+=", "\[Delta]f"}], ";"}], "}"}]}], "]"}]}], 
  "*)"}]], "Input",
 CellChangeTimes->{
  3.456746939*^9, {3.5372003913424997`*^9, 3.5372003937643747`*^9}, {
   3.5378885919057503`*^9, 3.5378885935932503`*^9}, {3.537901719773398*^9, 
   3.53790172663741*^9}, 3.5379017983975363`*^9, 3.5385235615778646`*^9, 
   3.538524008413422*^9, {3.540084146036523*^9, 3.5400841470661955`*^9}, 
   3.5404312983500743`*^9, 3.5404424331582637`*^9, {3.5406998704760923`*^9, 
   3.540699870788086*^9}, {3.5407000208571997`*^9, 3.5407000214031897`*^9}, {
   3.5407003545674286`*^9, 3.540700373396387*^9}, {3.5407005292071896`*^9, 
   3.54070054245142*^9}, {3.540701011319009*^9, 3.5407010167789383`*^9}, {
   3.540701078148552*^9, 3.54070107905334*^9}, {3.540702894168069*^9, 
   3.54070290344995*^9}, 3.540704212332471*^9, 3.54070429143855*^9, {
   3.5407339187505383`*^9, 3.5407339258319435`*^9}, {3.5407346581998324`*^9, 
   3.5407346626740885`*^9}, {3.5407390969557147`*^9, 3.540739098918827*^9}, 
   3.540740429510459*^9, {3.5410306331138387`*^9, 3.5410306802726383`*^9}, {
   3.5410308040274386`*^9, 3.541030814635439*^9}, {3.5412704895519075`*^9, 
   3.5412704916584883`*^9}}],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"rFrep", "=", "0.0"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"tt", "=", 
    RowBox[{"AbsoluteTime", "[", "]"}]}], ";"}], " ", 
  RowBox[{"(*", 
   RowBox[{
   "distribute", " ", "NIntegrate", " ", "among", " ", "CPU", " ", "cores"}], 
   "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"specList", "=", 
   RowBox[{"{", "}"}]}], ";", 
  RowBox[{
  "SetSharedVariable", "[", "specList", "]"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"ParallelDo", "[", 
   RowBox[{
    RowBox[{"AppendTo", "[", 
     RowBox[{"specList", ",", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{
         "fRepList", "\[LeftDoubleBracket]", "kk", "\[RightDoubleBracket]"}], 
         ")"}], ",", 
        FractionBox[
         RowBox[{"Chop", "[", 
          RowBox[{"vIntSpec", "[", 
           RowBox[{
            RowBox[{
            "fRepList", "\[LeftDoubleBracket]", "kk", 
             "\[RightDoubleBracket]"}], ",", 
            RowBox[{
            "fOffsetList", "\[LeftDoubleBracket]", "kk", 
             "\[RightDoubleBracket]"}]}], "]"}], "]"}], "12"]}], "}"}]}], 
     "]"}], ",", 
    RowBox[{"{", 
     RowBox[{"kk", ",", 
      RowBox[{"Length", "[", "fRepList", "]"}]}], "}"}]}], "]"}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"ttt", "=", 
   RowBox[{
    RowBox[{"AbsoluteTime", "[", "]"}], "-", "tt"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{"Print", "[", 
  RowBox[{"\"\<Execution time = \>\"", ",", "ttt"}], "]"}]}], "Input",
 CellChangeTimes->{{3.540700569251876*^9, 3.540700631292281*^9}, {
   3.540704216154398*^9, 3.540704252064907*^9}, {3.5407042969296446`*^9, 
   3.5407042987392097`*^9}, 3.5407055725228095`*^9, {3.54073392927214*^9, 
   3.540734037008302*^9}, {3.541275778337726*^9, 3.541275794685014*^9}, {
   3.541276893218398*^9, 3.541276895647437*^9}, {3.541296195763912*^9, 
   3.541296442665784*^9}, {3.5412965326870537`*^9, 3.5412965488659134`*^9}, {
   3.541298036226572*^9, 3.5412980663070354`*^9}, {3.5412999438710012`*^9, 
   3.5413000247309847`*^9}, 3.5413009553674545`*^9, {3.541302060074316*^9, 
   3.54130207288249*^9}, {3.541302810233549*^9, 3.5413028515127363`*^9}, {
   3.541304182880085*^9, 3.54130419865219*^9}, {3.542556980744997*^9, 
   3.542556983272213*^9}, {3.5425636106215625`*^9, 3.542563611370367*^9}, {
   3.54922355984375*^9, 3.5492235624853516`*^9}}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"Execution time = \"\>", "\[InvisibleSpace]", 
   "2010.491211`10.754847172390987"}],
  SequenceForm["Execution time = ", 2010.491211`10.754847172390987],
  Editable->False]], "Print",
 CellChangeTimes->{3.549217118404297*^9, 3.549221008463867*^9, 
  3.549223402029297*^9, 3.549225609551758*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{
    RowBox[{"totalSpec", "=", 
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{
         RowBox[{
         "fRepList", "\[LeftDoubleBracket]", "t", "\[RightDoubleBracket]"}], 
         ",", 
         RowBox[{"Chop", "[", 
          RowBox[{
          "specList", "\[LeftDoubleBracket]", "t", "\[RightDoubleBracket]"}], 
          "]"}]}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"t", ",", "1", ",", 
         RowBox[{"Length", "[", "specList", "]"}]}], "}"}]}], "]"}]}], ";", 
    "\n", 
    RowBox[{"fRepList", "=", 
     RowBox[{"fRepList", "-", 
      RowBox[{"Min", "[", "fRepList", "]"}]}]}], ";", "\[IndentingNewLine]", 
    RowBox[{"plotSpec", "=", 
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{
         RowBox[{
         "fRepList", "\[LeftDoubleBracket]", "t", "\[RightDoubleBracket]"}], 
         ",", 
         RowBox[{"Chop", "[", 
          RowBox[{
          "specList", "\[LeftDoubleBracket]", "t", "\[RightDoubleBracket]"}], 
          "]"}]}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"t", ",", "1", ",", 
         RowBox[{"Length", "[", "specList", "]"}]}], "}"}]}], "]"}]}], ";"}], 
   "*)"}], "\[IndentingNewLine]", 
  RowBox[{"ListPlot", "[", 
   RowBox[{
    RowBox[{"Sort", "[", "specList", "]"}], ",", 
    RowBox[{"PlotRange", "\[Rule]", "All"}], ",", 
    RowBox[{"Joined", "\[Rule]", "True"}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.5409271816468678`*^9, 3.5409272231595325`*^9}, 
   3.540959220519678*^9, {3.541299917754927*^9, 3.5412999284104104`*^9}, {
   3.5413019390284843`*^9, 3.5413019427882533`*^9}, {3.5413033228068624`*^9, 
   3.5413033265042048`*^9}, {3.541303499156045*^9, 3.5413034998112698`*^9}, {
   3.5492213639384766`*^9, 3.549221367432617*^9}, {3.5492235446572266`*^9, 
   3.549223547359375*^9}}],

Cell[BoxData[
 GraphicsBox[
  {Hue[0.67, 0.6, 0.6], LineBox[CompressedData["
1:eJwl2Hk8VtseBvBNMylpOJqEQ4PKoW5KMlUIlUjSZApNNKgoFW8SpWgyTykl
hRDdkKFEhqKJSEolJVFoOEdHvbfnuX/1+X6e1lp7r73Wb62XktM2KxdJQRA2
SggC/i1VFAkuVWmGTt67Zj19eVG/VFGQgKPz6q8UJl7+bRH9rltGMbcwTf+o
oSAJlzu2T9YoztTXafKluxT39YZVXUf7PnC95eGX9oNv6guCiNYU1ruPGHnr
9/8X067m041O77iNvC+coXNDsPzn1m8L/WAf5UynItsi5HT03SwLpSt5aE8H
6R8+Xed+E+P3Z38Vdz888GFO71k1Rc8gLw/PPwC21PB75rIhHzlt+GOvsUls
EdoPhBMSYx+lhxVj/EHwsyS5McFVd/W9Cn3pC4umtFd030cuBXd/XNlma/cY
z0e/sgiUnSnxFP3Tj+pTLcKz6tC/NBz54Eep8YY6tB8Mj95/8HarJ/6/L+05
c4auSp9a5DLw1uGHxFcCHyOnLY2W/3tU9wHmdwjfz0IyTNqvEuPTx3/ZqXUp
lmF8OrJrcaKsTSnyoXDYbK/JGy6VoL0sXJbdG7NZuxT903uVK3PtZ5Rh/GGw
ywPn24veVmL+6FdjAhd+dH+IXA5+qPf0hEVrDXI6+GWDV2BcA/LhcM28znH/
FjRhfHrJpl7t7N2vkI/g9x7g+nKq9gvktNqWiR1BmZgvMd12T67leQPndyRc
qvhjaPGBajz/KHriwskRmyuQ0/LNkltbi/E+Yjo77bLJXKkS5H/A196Pvrqr
nt9XHl4srdIgPQMW0TYKQ85W9WA9CqP5/10uecf94vqk7W68nzWg8jbGHwMP
9F3hfM6ROe0TuCXo+5LbGJ+eE+Oe69vJfCwc4vbAKHLnHbQfx/69NVWjn/L7
0OsW5VVf1cH7CONhoWpfyPRBeF8RHbm4+0zAFMyHoAD3y3sa0OpYi5zu+eom
1fmZ62sC3JBueW1Y2DP0T9dmuThcV3qOXJG5ifP7haefoz2d6VJV7DOB30+J
/XUqrlmZyPVL98sVuzQdrv29H8S05FCx4d3UGuTK7P/sxwcKu2vx/vR26yHR
4ZL16P9PeNoc+fL9FxvRvwocMWFf0LCs12hPu+U5Lbce/g65Kr9n3Lt37kUf
kNOnHHx6pOw/Ip/I+qKUW7D0+kfsT1re73nmsrKPGJ++a/7p0iWXNow/ifNX
UK5WlfgO8z8Zfqsv1m/6/hrzQ+upR6V05GM+hCnwVf3jSod2PEV7OtbI2D44
A/tRTM8JOP1muw32g0gNDp/dnrpGvhrPNxUekRvsMs+rCjm96JvpyYuNzKfB
cz+4h8urPEJOf2ysWF04oQb5dFixss6kum89cvqtbfvEJpXnmH96UXGucu17
fC+ROnxlz8jykRPq8Px/wdrtlyek23L/0KoSQ/qZFt5HezpS11iyYSj2u0gD
Xiszo1J60h28Hy0nVz/ojQbqta8mrOA1ZtGiZaiXYrrFf/fxlCcFaD+D+9+h
/4GOsgLktIH6l8wns1i/Z8LxG1P+3inH/fMfuL+2cUzzSOwHEW11qjOkrqIE
58ssuG/DrL80xSUYnzZONbmhfeYO5kcLTh0uMchzJ/cXLZO+sWfcqCL0P5v7
uaBzeUp2PnL6vIW3SX4eziNhDmwUr5h84v1NjE9nfAn5nH6e5wNtoXUobc9A
vr826/O1jgHZpsXofy5cOd9q0uT8UrSnjf7R0ep5WIn2dJTVrh9Hah6ivQ4s
ZWh5fk476qUwD1bWdk23U+H3pY2Vbmk2zMZ+EXR5ntnUyMXH4TwQ0ZmNAzxj
crj/9eDdo9aYn9jD/ukcr/lTUodzfenD/rOaIk4Zcv3RugefyJYsY/00gEeO
SrpfEsr1S0etkA7KcHqC9oZwqI5BSnA16omIble4ZXTR8SXaz4f/WlJqdO9c
M8anaw9Lf9i4oxX5Aq6foEqn4QuxH33pZ1NMHHv029D/Qvhwy9nB39ux30V0
xbv04xlLm5Eb0UW/5kbNwXkhomeNv28ZuIT1yxgeJznEITyP9Yl2ct5pY5zH
89MEvnUj7l6pK+sjXZXh7r8u+SnyRbD64INbjB+yPtJRBpuior++RG4KO6ds
3r1Kqhnt6U0/B9/Musb6Ycb6s2LS9YiUVoz/f38ttZ++HRbTkqeKh7zVe4/c
HH4xd+zcm7dbsP/ondJ1juKH7H8x11eaYpTte9RDMa2zuntQgBPffwnXS5N9
4b0wnJfCUvhggp75vst8fnphzNmeR66oV4IFnLxynZ+fNeqJiF42v6j3W+oT
9E+PmqwQOGIa82VwWE7LZ3VtzJdgCX+3PrYnYRzrHz13kJe/yhPeb2gPK6O3
+ZtQb0RW8OOf9r1dA3l+0yPS5ZVmmnD9LodLmuOUH7/l+rKm+/XMNB/C+w2d
tEXxQKwh7xcr4OOzG7UPy1ehf1r2qU3O8I+4jwg28H7pG8OkqirQns5Ijjjz
fUk52q+Ec+c/tl7gfBfj0+Xa4oCWeN5/bLm+t12eFrsS9c6XHjP/q/Pfv+8H
v78P7Vnz/k33rWLkq+CJy0OOPTvD+8NqePXU5Yfvr2JOB9kFGDXkwWJaSzdO
w/s16pVoDfdrsfwe13LUC2EtfOGZ86Ysmwq8H327TqTm/if37zr4+pgZRyTW
P0JOf/WcV/9HFPe/HfvfWeUU+ID1nX7UJj5S/OIx+rfnefTFfl/nBs4vXeRv
4S97/iFyB3i0g7ymngzusyL6uDD5ZrddGXJHntemFg23lXl/ood4SXgst2P9
c4J/ao1Jl/zK+kvXWKqMXN/MfD084Eeh25UW1HMRveRS0VTHKt5P6U8OH8KH
/+L9zZn7V01tobom7rOCC7/f9ct7m9WZ037T9Jp3yPD7uvL/V/pYbbLm96U9
tnSFKbjyfrWB9a9+jqnqxFuYPzpKdbZ1yCDW/408D0TKodOW3kB72jP5avmB
zmzkm+D0qy0eBbFZyOm5T4+mnTeDhc2wQ8dVzcakLPRPq21v6TlgnY33o522
a+k3R2Yj3wLLP1ofWj8GFtxYL7737BOMsnB/oV+NDlg4Qfka+neHd5k3je3+
zN9XtJ2GhPySPrCYnhJyufyNfgbyrXDIizAh9UI62m9jfa28tVhyNSyi537x
/JJ5kvl2+KH+zsUqmmxPb/HKSQsbmIl8BxxzebzM9u8cn/5hbzbeXgfvK3jw
vJdNmPUmiu9Hn/b4e+z1P/6L9jt5fps+zxWnwSK6cKjHm5LXOXh+enqZ+tLK
SQVov4vzrbN4S+IS/v7bDTduy1sWPIb7i9aYaP26tgzrUfDkfDmNaO86VI6c
3tUWWT+7jb/PvFj/3b29DNZy/9APzPLveBrwfNgD/0c6oMUtnPcnetieQqWv
UrgPCnvhuMZtEY6uPJ/ox199ZpvLfkLuzfWclKFQpPsVOW25aqrTT+0fmJ99
cELk2JVznX5h/dNOSflZfSJ/It8PC1/qnz3P+Ac5bX4qYJ3dqC/ID3B/mCQO
mGnXgf7pqhVhA5Md32N8H86nkkf/IReb0Z6Oc8lPkY7n8/tyvwoNNat/8vyg
x63/npieylwEd/y6IhOZ9gbzRxck2stevMT+D8LdazRyh8zuQE57ezcX5zh3
Iffj/HkcSj37326sX9pCOvNT8LQu3N8OwaVTX2SFHfiI+acTVKtj06Ra8P3p
hr37b1yX5vP5w26bzio7jMN5Jhzm99o7NfCcM88neqzZ0U4NQ9yXhQBYodUg
JW5zNZ6PPpr0yu/EYn7/QHhj+w7Zwpx7GJ+WVO9naTcS9yHhCNt3l7ltPPsA
/dOqeSvDX4p5/hzl/t5gEr+nkPcj+vSDE7r77arwfkFwV8jPNSefo16K6LbH
2xrPjeD9mb4SHnTwQN8izM8xnjdOBp8N63OR0zLuradDNbk/jsNp45J1lLxR
D8S0nFuBSj8v7EdRMJxzXnNAZwf3bwicv1lL+2RRGnLabnTp0jSzFDzfCa7X
kC4jxxb+/YbukPlk9DoqGevrJGxqs/bDxcpkzB9tpn/ST3dYEtqfgg3GznTw
eJGI/un7mydq+/qfx/PTg765Bp9vO4v3O836N07HYV17Avo/w/PReamEthws
ouVrI/1lQ+Px/KHcr9fc49s2xGF8umR8iPndkFjkYXCY07HtiWYxaE9XvupV
/ssnGvNDP9aarrn+XhTycK532422WzVhIQKuXRFd0TsnCv3TVgsCR/Wri0T/
kaxnkVlJEYGwiM4cIRnZNpN5FNdHtEqzpypz+m7cu3PJAyPRfzS8aFV23bFJ
keifjnR7+yRiG9vHwEf+1Ja/eTIS80MvNdxo1n9KFPJY1usJ6krK4Xw+enWI
5DTdO3z+ONhFa8q2Oslo5HSvrOnwTY1sH8/zLkGhz5qyaLw/nfSysXdCMuZL
OMv1sv2SU6goDs9PR4fdUR+Tz/lP4Pngm6HbxzkBOd0Ssd/vpya/L22Ts/uA
ac4FjH8Obrc6Wrm97RLan4e/ra1pUp+XivZ09o5PHVILMpAncj1U3HFrleb6
pc3U2p4d2Z2B9XWB3/fzkgWKX68ip/0UQl8/jk/D81/k8+ofqQ2y5fqmWyd2
9erLXcXz0SUrAtb5y/L8SoLFC103FJ7j+UjrB4ZP3FOE81x0CY7/1/9PGcVc
PF8yXH2rLrFqaxHqJ623MVixbjJ+b4rp1p6eGCtn3s8uc722Gmls2MH6cQXW
LV86p3Mo7msiOs0yd/DxL/eQp7C+FxeustvK9rRe39Bij1TcV4RU9vfy++eo
c7wf0j8lD8f11uO+IqRxf1WH6b2J4u9rWuJA10pVVfy9U0y/fjLZyeePPORX
uR5zfP5e5wYL6ZzPUnc/Vfd8jE+HV4nnBI7n31/pYFOtUadmc/wMng/eI8YV
9sV5Kqa/FRislKmrxvxkwhf/UdWo2sn7+zXYdu0RNf8Y/v2SVvQ0a/GL4u+P
LM7vc9VZOWeeYP3THpJrTA9uQ70VsnleekrHqjeV6/8PzFNfbw==
    "]]},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->True,
  PlotRangeClipping->True]], "Output",
 CellChangeTimes->{3.549225610213867*^9}]
}, Open  ]],

Cell[BoxData[""], "Input",
 CellChangeTimes->{{3.5492245799990234`*^9, 3.5492245800751953`*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Export", "[", 
  RowBox[{"oFile", ",", 
   RowBox[{"Sort", "[", "specList", "]"}], ",", "\"\<Table\>\""}], 
  "]"}]], "Input",
 CellChangeTimes->{{3.541303835300945*^9, 3.5413038495754023`*^9}, {
  3.549221369826172*^9, 3.5492213706152344`*^9}, {3.5492235499697266`*^9, 
  3.5492235506347656`*^9}}],

Cell[BoxData["\<\"C:\\\\Users\\\\work\\\\Data\\\\output\\\\Rb87-D1.txt\"\>"], \
"Output",
 CellChangeTimes->{3.5492171192148438`*^9, 3.5492210092226562`*^9, 
  3.5492234027978516`*^9, 3.5492256103154297`*^9}]
}, Open  ]]
}, Open  ]]
},
WindowSize->{1284, 785},
WindowMargins->{{Automatic, -407}, {-155, Automatic}},
CellContext->Notebook,
FrontEndVersion->"8.0 for Microsoft Windows (32-bit) (November 7, 2010)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[545, 20, 460, 13, 72, "Input"],
Cell[CellGroupData[{
Cell[1030, 37, 76, 1, 49, "Subtitle"],
Cell[1109, 40, 4015, 91, 292, "Input"],
Cell[CellGroupData[{
Cell[5149, 135, 37, 0, 28, "Subsubtitle"],
Cell[5189, 137, 148, 5, 31, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[5374, 147, 40, 0, 28, "Subsubtitle"],
Cell[5417, 149, 223, 8, 31, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[5677, 162, 107, 1, 28, "Subsubtitle"],
Cell[5787, 165, 330, 10, 31, "Input"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[6166, 181, 75, 1, 49, "Subtitle"],
Cell[6244, 184, 102, 1, 49, "Subtitle"],
Cell[6349, 187, 149, 2, 28, "Subsubtitle"],
Cell[6501, 191, 325, 7, 47, "Input"],
Cell[CellGroupData[{
Cell[6851, 202, 105, 1, 28, "Subsubtitle"],
Cell[6959, 205, 408, 10, 47, "Input"],
Cell[CellGroupData[{
Cell[7392, 219, 105, 1, 31, "Input"],
Cell[7500, 222, 161, 2, 30, "Output"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[7710, 230, 101, 1, 28, "Subsubtitle"],
Cell[7814, 233, 219, 6, 31, "Input"],
Cell[8036, 241, 167, 5, 31, "Input"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[8252, 252, 106, 1, 49, "Subtitle"],
Cell[8361, 255, 222, 3, 28, "Subsubtitle"],
Cell[8586, 260, 306, 7, 47, "Input"],
Cell[CellGroupData[{
Cell[8917, 271, 199, 3, 28, "Subsubtitle"],
Cell[9119, 276, 287, 7, 47, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[9443, 288, 109, 1, 28, "Subsubtitle"],
Cell[9555, 291, 301, 8, 52, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[9893, 304, 389, 5, 28, "Subsubtitle"],
Cell[10285, 311, 925, 17, 31, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[11247, 333, 413, 6, 28, "Subsubtitle"],
Cell[11663, 341, 554, 9, 31, "Input"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[12266, 356, 180, 4, 49, "Subtitle"],
Cell[12449, 362, 249, 3, 28, "Subsubtitle"],
Cell[12701, 367, 330, 7, 47, "Input"],
Cell[CellGroupData[{
Cell[13056, 378, 199, 3, 28, "Subsubtitle"],
Cell[13258, 383, 403, 8, 47, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[13698, 396, 115, 1, 28, "Subsubtitle"],
Cell[13816, 399, 347, 9, 52, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[14200, 413, 252, 3, 28, "Subsubtitle"],
Cell[14455, 418, 630, 16, 47, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[15122, 439, 315, 4, 28, "Subsubtitle"],
Cell[15440, 445, 1230, 22, 52, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[16707, 472, 350, 5, 28, "Subsubtitle"],
Cell[17060, 479, 1058, 19, 31, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[18155, 503, 371, 5, 28, "Subsubtitle"],
Cell[18529, 510, 1053, 19, 31, "Input"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[19631, 535, 155, 2, 49, "Subtitle"],
Cell[19789, 539, 240, 3, 28, "Subsubtitle"],
Cell[20032, 544, 398, 8, 47, "Input"],
Cell[CellGroupData[{
Cell[20455, 556, 199, 3, 28, "Subsubtitle"],
Cell[20657, 561, 473, 9, 47, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[21167, 575, 104, 1, 28, "Subsubtitle"],
Cell[21274, 578, 347, 9, 52, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[21658, 592, 293, 4, 28, "Subsubtitle"],
Cell[21954, 598, 764, 16, 47, "Input"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[22767, 620, 139, 1, 49, "Subtitle"],
Cell[22909, 623, 175, 2, 28, "Subsubtitle"],
Cell[23087, 627, 442, 15, 31, "Input"],
Cell[23532, 644, 1602, 47, 56, "Input"],
Cell[25137, 693, 2201, 51, 212, "Input"],
Cell[CellGroupData[{
Cell[27363, 748, 57, 0, 28, "Subsubtitle"],
Cell[27423, 750, 1384, 45, 76, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[28844, 800, 97, 1, 28, "Subsubtitle"],
Cell[28944, 803, 1254, 33, 99, "Input"],
Cell[30201, 838, 1184, 32, 99, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[31422, 875, 115, 1, 28, "Subsubtitle"],
Cell[31540, 878, 2134, 54, 275, "Input"],
Cell[33677, 934, 2166, 55, 275, "Input"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[35892, 995, 89, 1, 49, "Subtitle"],
Cell[35984, 998, 816, 19, 74, "Input"],
Cell[CellGroupData[{
Cell[36825, 1021, 280, 7, 53, "Input"],
Cell[37108, 1030, 165, 2, 30, "Output"]
}, Open  ]],
Cell[37288, 1035, 222, 7, 50, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[37547, 1047, 73, 1, 49, "Subtitle"],
Cell[37623, 1050, 39, 0, 28, "Subsubtitle"],
Cell[37665, 1052, 794, 15, 52, "Input"],
Cell[CellGroupData[{
Cell[38484, 1071, 102, 1, 28, "Subsubtitle"],
Cell[CellGroupData[{
Cell[38611, 1076, 1491, 24, 52, "Input"],
Cell[40105, 1102, 159, 2, 30, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[40301, 1109, 206, 5, 47, "Input"],
Cell[40510, 1116, 158, 2, 30, "Output"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[40717, 1124, 90, 1, 28, "Subsubtitle"],
Cell[40810, 1127, 558, 16, 87, "Input"],
Cell[41371, 1145, 833, 25, 87, "Input"],
Cell[CellGroupData[{
Cell[42229, 1174, 1267, 27, 31, "Input"],
Cell[43499, 1203, 407, 8, 30, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[43943, 1216, 1321, 28, 31, "Input"],
Cell[45267, 1246, 336, 7, 30, "Output"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[45652, 1259, 102, 1, 28, "Subsubtitle"],
Cell[45757, 1262, 391, 10, 31, "Input"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[46197, 1278, 112, 1, 49, "Subtitle"],
Cell[46312, 1281, 6012, 148, 283, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[52361, 1434, 185, 3, 49, "Subtitle"],
Cell[52549, 1439, 2031, 40, 172, "Input"],
Cell[CellGroupData[{
Cell[54605, 1483, 2413, 60, 148, "Input"],
Cell[57021, 1545, 343, 7, 23, "Print"]
}, Open  ]],
Cell[CellGroupData[{
Cell[57401, 1557, 1885, 49, 92, "Input"],
Cell[59289, 1608, 5299, 90, 242, "Output"]
}, Open  ]],
Cell[64603, 1701, 96, 1, 31, InheritFromParent],
Cell[CellGroupData[{
Cell[64724, 1706, 322, 7, 31, "Input"],
Cell[65049, 1715, 208, 3, 30, "Output"]
}, Open  ]]
}, Open  ]]
}
]
*)

(* End of internal cache information *)
